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ABSTRACT 

SN 1997ef has been recognized as a peculiar supernova from its light curve and 
spectral properties. The object was classified as a Type Ic supernova (SN Ic) because 
its spectra were dominated by broad absorption lines of oxygen and iron, lacking any 
clear signs of hydrogen or helium line features. The light curve is very different from 
that of previously known SNe Ic, showing a very broad peak and a slow tail. The 
strikingly broad line features in the spectra of SN 1997ef, which were also seen in the 
hypernova SN 1998bw, suggest the interesting possibility that SN 1997ef may also be 
a hypernova. 

The light curve and spectra of SN 1997ef were modeled first with a standard SN Ic 
model assuming an ordinary kinetic energy of explosion £k = 10 51 erg. The explosion 
of a CO star of mass Mco ~ 6M gives a reasonably good fit to the light curve but 
clearly fails to reproduce the broad spectral features. Then, models with larger masses 
and energies were explored. Both the light curve and the spectra of SN 1997ef are much 
better reproduced by a C+O star model with = 8 x 10 51 erg and Mqq = WMq. 
Therefore, we conclude that SN 1997ef is very likely a hypernova on the basis of its 
kinetic energy of explosion. Finally, implications for the deviation from spherical 
symmetry are discussed in an effort to improve the fits to the observations. 

Subject headings: gamma rays: bursts — radiation transfer: light curves — radiation 
transfer: spectra — stars: massive — stars: supernovae — stars: supernovae: 
individual (SN 1997ef) 
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1. INTRODUCTION 

The supernova 1997ef (SN 1997ef) was discovered on November 25, 1997 at an R magnitude 
of 16.7 near the spiral galaxy UGC4107 (Sano 1997). The first spectrum was taken on November 
26 (Garnavich et al. 1997a). Subsequently, photometric and spectroscopic follow-ups have 
provided high quality optical light curves and spectra (Garnavich et al. 1997a,b,c; Hu et al. 1997; 
Filippenko et al. 1997; Wang et al. 1998). As seen in Figure H the spectra of SN 1997ef are 
dominated by broad oxygen and iron lines but do not show any clear feature of hydrogen or helium 
(Garnavich et al. 1997b; Filippenko et al. 1997), showing the overall similarity to other Type Ic 
supernovae (SNe Ic) 19941 and 1998bw. This led us to classify SN 1997ef as a SN Ic. 

In Figure |2| the visual light curve of SN 1997ef (Garnavich et al. 1997b, c) is compared with 
those of the SN Ic 1998bw (Galama et al. 1998) and the ordinary SN Ic 19941 (Richmond et al. 
1996a, b). Despite the spectral similarity, the light curve of SN 1997ef is quite different from those 
of SN 1998bw and SN 19941. It has quite a flat peak, much broader than those of the other SNe Ic. 
Besides, the tail of the light curve of SN 1997ef starts late and the rate of its decline is much slower 
than in other SNe Ic. It is also true that the light curves are rather diverse, even in this limited 
number of samples, implying a range of energies and/or progenitor masses of SN Ic explosions. 

The most striking and peculiar characteristic of SN 1997ef is the breadth of its line features. 
Such broad spectral features were later recognized to be a distinguishing property of the spectra in 
SN 1998bw (Fig. |). SN 1998bw was discovered within the error box of GRB980425 determined 
by the BeppoSAX satellite, only 0.9 days after the date of the gamma-ray burst (GRB), and 
therefore probably related to this GRB (Galama et al. 1998). The very broad spectral features 
and the light curve shape have led to the conclusion that SN 1998bw had an extremely large 
kinetic energy of explosion, Ek ~ 3 x 10 52 ergs (Iwamoto et al. 1998; Woosley, Eastman, & 
Schmidt 1999). This was one order of magnitude larger than the energy of typical supernovae, 
thus SN 1998bw was termed a "hypernova" (Iwamoto et al. 1998). 

The spectral similarities between SN 1997ef and SN 1998bw suggest the interesting possibility 
that SN 1997ef may also be a hypernova. In fact, a possible connection with a GRB has been 
suggested for SN 1997ef: GRB971115 appears to be compatible with the supernova in the position 
and the time of occurrence (Wang Sz Wheeler 1998). Since the statistical significance for this case 
is much weaker than for the case of SN 1998bw and GRB980425, it is difficult to confirm the 
physical association between SN 1997ef and GRB971115. However, it is possible at least to clarify 
whether SN 1997ef is a hypernova or not by estimating the kinetic energy of explosion through 
modeling of light curves and spectra as in the case of SN 1998bw (Iwamoto et al. 1998; Mazzali 
1999). This is exactly the primary purpose of this paper. 

We constructed supernova progenitor models and performed detailed hydrodynamics and 
radiation transfer calculations to obtain light curves and spectra for the explosion models. The 
results were compared with observations of SN 1997ef, to derive its explosion energy and the 
ejecta mass, and thus to determine whether SN 1997ef was an ordinary SN Ic or a hypernova. 
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Since the light curves of the other SNe Ic 19941 and 1998bw were successfully reproduced by the 
collapse-induced explosion of C+O stars (Nomoto et al. 1994; Iwamoto et al. 1994; Iwamoto et 
al. 1998), we adopted C+O stars as progenitor models for SN 1997ef as well. 

This paper consists of six sections including this introduction. Section 2 describes the ordinary 
SN Ic model and the hypernova model. The method and results of our light curve calculations 
are presented in §3. The synthetic spectra are compared with the observations in §4. Section 5 is 
devoted to discussion on various issues including the SN-GRB connection and possible progenitor 
scenarios. Finally, our conclusions are summarized in §6. 



2. EXPLOSION MODELS FOR SUPERNOVAE AND HYPERNOVAE 

We construct hydrodynamical models of an ordinary SN Ic and a hypernova as follows. 

(1) In the ordinary SN Ic model(model CO60), a C+O star with a mass Mqo = 6.OM 
(which is the core of a 25 M main-sequence star) explodes with kinetic energy of explosion 
E K = 1.0 x 10 51 ergs and ejecta mass M cj = M C o ~ M rcm = 4.6M . Here M rem (= 1.4 M ) 
denotes the mass of the compact star remnant (either a neutron star or a black hole). 

(2) In the hypernova model (CO100), a C+O star of Mqo = 1OM is constructed from 
the 10 Mq He star (which has a 8 M C+O core) by removing the outermost 2 M of He 
layer and extending the C+O layer up to 1OM . This model corresponds to 30 - 35 M on the 
main-sequence. This progenitor goes off with Ek = 8.0 x 10 51 ergs and M c j = 7.6M , i.e., M rem = 
2.4 M Q . 

The hydrodynamics at early phases was calculated by using a Lagrangian PPM code (Colella 
& Woodward 1984) with a simple nuclear reaction network including 13 alpha elements (Miiller 
1986). Detailed post-processing calculations were carried out with a larger-size nuclear reaction 
network including 240 isotopes (Hix & Thielemann 1996). The explosion is triggered by depositing 
thermal energy in a couple of zones just below the mass cut so that the final kinetic energy 
becomes the required value. The position of the mass cut is adjusted for the ejected mass of 56 Ni 
to be M( 56 Ni) = O.15M . 

The compact remnant in CO60 is likely a neutron star because M rcm = 1.4 M , while it 
may be a black hole in CO100 because M rcm (= 2.4 M ) may well exceed the maximum mass 
of a stable neutron star. The above values of M rem are determined so that M ( 56 Ni)= O.15M is 
ejected to reproduce the maximum brightness of SN 1997ef by the radioactive decay heating of 
56 Ni and 56 Co. 

These model parameters are summarized in Table 1. They can be constrained by comparing 
the calculated light curves and synthetic spectra with observations. The parameters of models 
C021 for SN 19941 (Nomoto et al. 1994; Iwamoto et al. 1994) and C0138 for SN 1998bw 
(Iwamoto et al. 1998) are also given in Table 1. We constructed the progenitor model by attaching 
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a thin hydrostatic and in-thermal-equilibrium C+0 envelope to the C+O core of the presupernova 
model (Nomoto & Hashimoto 1988; Hashimoto 1995). 

The expansion soon becomes homologous so that v oc r. The solid lines in Figure [| show the 
density distributions in the velocity space for CO60 and CO100 at t = 16 days. The expansion 
velocities are clearly higher in CO100 than in CO60. Figures |] and || show the composition 
structure of models CO60 and CO100, respectively, against the expansion velocity and the 
Lagrangian mass coordinate of the progenitor. In CO100, the Fe and Si-rich layers expand much 
faster than in CO60. The total amount of nucleosynthesis products are summarized in Table 2. 



3. LIGHT CURVE MODELS 

3.1. Radiation Hydrodynamics Code 

The light curve calculation is started, when the ejecta have reached the homologous expansion 
phase, with a one-dimensional spherically symmetric radiation transfer code (Iwamoto 1997). 

The code solves the multi-frequency radiative transfer equation for the specific intensity l v in 
the comoving frame, including all terms up to the first order in v/c (Mihalas & Mihalas 1984), 
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where k v and o v are the absorptive and scattering opacities, respectively, B u is the Planck function 
and ti is the cosine of the angle made by the radial direction and the direction of the ray. This 
equation is solved numerically using the Feautrier method with an approximate Lambda operator 
similar to the one described by Hauschildt (1992). 

To determine the gas temperatures, Equation (1) is solved simultaneously with the energy 
equation and the first two moment equations of Equation (1). The energy equation of the radiation 
plus gas is written as 
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while the radiation energy and momentum equations are 
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respectively, where e is the thermal energy of ions and electrons per unit mass, and E, F and / 
are the radiation energy density, flux, and the Eddington factor defined as follows. 
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Partial derivatives with respect to t in Equations (l)-(4) are all Lagrangian time derivatives. 
The absorptive and scattering parts of the opacity are given as 

K v = e(K b _ b + K b _f) + Kf_f , (8) 

and 

OV = (1 - e)(«b-b + K b-f) + ^eO"T, (9) 

where Kb-b 5 K b-f and Kf_f are the bound-bound, bound-free, and free- free opacities, 
respectively; n e is the number density of free electrons, and ctt is the Thomson scattering cross 
section. In the moment equations, the energy mean (ke) and the Planck mean opacities (rep) 
include only the absorptive part, while the flux mean opacity (xf) is the total opacity. 

For bound-bound transitions, energy levels and transition probabilities are taken from the 
compilation by Kurucz (1991). For bound-free data, we use the analytic fitting formula to 
the photoionization cross sections given by Verner & Yakovlev (1995). Local Thermodynamic 
Equilibrium (LTE) is assumed to determine the ionization balance and the level populations of 
each ion. However, the non-LTE effect is approximately taken into account by assuming that 
the value of the absorptive fraction e is a constant less than unity in Equations (8) and (9). 
Experiments of spectral syntheses have shown that a value e = 0.1 is a reasonable choice for this 
fraction (Baron et al. 1996). 

It has been argued by several authors that the expansion would increase the chance of 
interactions between radiation and matter through line transitions, and thus the mean opacities 
become larger than in static medium, especially at relatively early phases (Karp, Lasher, & Chan 
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1977; Eastman & Pinto 1993; Blinnikov 1996). However, we neglected this expansion effect in 
evaluating the mean opacities, although this is one of the currently controversial issues. 

The energy deposition due to the radioactive decays is calculated with a one energy-group 
7-ray-transfer code (Iwamoto 1997). We assume an absorptive opacity k 7 = 0.03 and the complete 
trapping of positrons. The rest-frame flux is calculated from the comoving-frame intensities using 
the following transformation: 



dfi, (10) 
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For the calculation of the light curves of CO60 and CO100 discussed in the next subsection, 
we use about 200 radial mesh points to solve the moment Equations (2)-(4), while 800 frequency 
and 50 radial mesh points were used for the multi-frequency radiative transfer Equation (1). 



3.2. Light Curve Models 

In Figure |6] we compare the calculated V light curves for models CO60 and CO 100 with the 
observed V light curve of SN1997ef. We adopt a distance of 52.3 Mpc (a distance modulus of 
fi = 33.6 mag) as estimated from the recession velocity, 3,400 km s -1 (Garnavich et al. 1997a) and 
a Hubble constant Hq = 65 km s _1 Mpc -1 . We assume no color excess (E(B — V) = 0.00); this 
is justified by the fact that no signature of a narrow Na I D interstellar absorption line is visible 
in the spectra of SN 1997ef at any epochs (Garnavich 1997a). The light curve of SN 1997ef has a 
very broad maximum, which lasts for ~ 25 days. This is much broader than in both the ordinary 
SN Ic 19941 and the hypernova SN 1998bw. The light curve tail of SN 1997ef starts only ~ 40 
days after maximum, much later than in other SNe Ic. 

The light curve of SN 1997ef can be reproduced basically with various explosion models with 
different energies and masses. In general, the properties of the light curve are characterized by the 
decline rate in the tail and the peak width, r pca k. The peak width scales approximately as 

r pcak oc ^M 3 J 4 E^\ (12) 

where k denotes the optical opacity (Arnett 1996). This is the time-scale on which photon diffusion 
and hydro dynamical expansion become comparable. Since the model parameters of CO 100 and 
CO60 give similar Tp ea k, the light curves of the two models look similar: both have quite a broad 
peak and reproduce the light curve of SN1997ef reasonably well (Fig. ||). 

The light curve shape depends also on the distribution of 56 Ni, which is produced in the 
deepest layers of the ejecta. More extensive mixing of 56 Ni leads to an earlier rise of the light 
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curve. For SN 1997ef, the best fit is obtained when the 56 Ni is mixed almost uniformly to the 
surface for both models. Without such extensive mixing, the rise time to V = 16.5 mag would be 
~ 30 d for CO100, which is clearly too long to be compatible with the spectroscopic dating (see 
§4). 

Model CO60 has the same kinetic energy (£"k = 1 x 10 51 erg) as model C021, which was 
used for SN Ic 19941 (see Table 1 for the model parameters). Since the light curve of SN 1997ef 
is much slower than that of SN 19941, the ejecta mass of CO60 is ~ 5 times larger than that of 
C021. 

The ejecta mass of CO100 is a factor of ~ 2 larger than that of CO60, and it is only ~ 20% 
smaller than that of model C0138, which was used for SN 1998bw (Table 1). Thus the explosion 
energy of CO100 should be ~ 8 times larger than that of CO60 to reproduce the light curve of 
SN 1997ef. This explosion is very energetic, but still much weaker than the one in C0138. The 
smaller £"k for a comparable mass allows CO 100 to reproduce the light curve of SN 1997ef, which 
has a much broader peak than that of SN 1998bw. 

The light curve of SN 1997ef enters the tail around day 40. Since then, the observed V 
magnitude declines linearly with time at a rate of ~ 1.1 x 10~ 2 mag day -1 , which is slower than in 
other SNe Ic and is even close to the 56 Co decay rate 9.6 x 10~ 3 mag day -1 . Such a slow decline 
implies much more efficient 7-ray trapping in the ejecta of SN 1997ef than in SN 19941. The ejecta 
of both CO 100 and CO60 are fairly massive and are able to trap a large fraction of the 7-rays, so 
that the calculated light curves have slower tails compared with C021. 

However, the light curves for both models decline somewhat faster in the tail than the 
observations. A similar discrepancy has been noted for the Type lb supernovae (SNe lb) 1984L 
and 1985F (Swartz h, Wheeler 1991; Baron, Young, & Branch 1993). The late time light curve 
decline of these SNe lb is as slow as the 56 Co decay rate, so that the inferred value of M is 
significantly larger (and/or £"k is smaller) than those obtained by fitting the early light curve 
shape. Baron et al. (1993) suggested that the ejecta of these SNe lb must be highly energetic and 
as massive as ~ 50 M . In §5.1, we will suggest that such a discrepancy between the early- and 
late-time light curves might be an indication of asphericity in the ejecta of SN 1997ef and that it 
might be the case in those SNe lb as well. 

3.3. Photospheric Velocities 

As we have shown, light curve modeling provides direct constraints on Mqq and However, 
it is difficult to distinguish between the ordinary SN Ic and the hypernova model from the light 
curve shape alone, since models with different values of M c j and can reproduce similar light 
curves. However, these models are expected to show different evolutions of the photospheric 
velocity and the spectrum as will be discussed in the following sections. 
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The photospheric velocity scales roughly as t> p h oc M c j Ej[ , so that M e j and Ek can be 
constrained by v p h in a different way from by means of the light curve width. Figure |7| shows the 
evolution of the observed velocities of the Si II line measured in the spectra at the absorption core, 
and the velocities at the grey photosphere computed by the light curve code for models CO60 and 
CO100. The velocities of the Si II line are somewhat higher than that of the photosphere, reaching 
~ 20,000 km s _1 at the earliest time. 

In model CO60 the photosphere forms at velocities much smaller than those of the observed 
lines, while CO 100 gives photospheric velocities as high as the observed ones. It is clear, from this 
comparison, that the hyper-energetic model CO100 is preferable to the ordinary model CO60. 
The apparent discrepancy that still exists between the CO 100 and observations might be related 
to the morphology of the ejecta, i.e., its deviation from spherical symmetry, as was also suggested 
in the case of SN 1998bw (Hoflich, Wheeler, & Wang 1999; Iwamoto et al. 1998). This issue will 
be discussed in §5.1. 



4. SYNTHETIC SPECTRA 

To strengthen the arguments in §3.3, we compare the observed spectra with theoretical 
model spectra computed using our explosion models with a more sophisticated spectrum synthesis 
code(Mazzali & Lucy 1993; Lucy 1999; Mazzali 1999). With such a detailed spectrum synthesis, 
we can distinguish between different models more clearly, because the spectrum contains much 
more information than a single-band light curve. 

Around maximum light, the spectra of SN 1997ef show just a few very broad features, and 
are quite different from those of ordinary SNe Ib/c, but similar to SN 1998bw. However, at later 
epochs the spectra develop features that are easy to identify, such as the Ca II IR triplet at 
~ 8200A, the O I absorption at 7500 A, several Fe II features in the blue, and they look very 
similar to the spectrum of the ordinary SN Ic 19941. 

We computed synthetic spectra with a Monte Carlo spectrum synthesis code using the density 
structure and composition of the hydrodynamic models CO60 and CO100. The code is based on 
the pure scattering code described by Mazzali & Lucy (1993), but has been improved to include 
photon branching, so that the reprocessing of the radiation from the blue to the red is followed 
more accurately and efficiently (Lucy 1999; Mazzali 1999). 

We produced synthetic spectra for three epochs near maximum, of SN 1997ef: Nov 29, Dec 5, 
and Dec 17. These are early enough that the spectra are very sensitive to changes in the kinetic 
energy. As in the light curve comparison, we adopted a distance modulus of \x = 33.6 mag, and 
E{B — V) = 0.00. The model parameters, the computed temperatures and the magnitudes of the 
synthetic spectra for CO100 are listed in Table 3. 

In Figure || we show the synthetic spectra computed with the ordinary SN Ic model CO60. 
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The lines in the spectra computed with this model are always much narrower than the observations. 
This clearly indicates a lack of material at high velocity in model CO60, and suggests that the 
kinetic energy of this model is much too small. 

Synthetic spectra obtained with the hypernova model CO 100 for the same 3 epochs are 
shown in Figure [| The spectra show much broader lines, and are in good agreement with the 
observations. In particular, the blending of the Fe lines in the blue, giving rise to broad absorption 
troughs, is well reproduced, and so is the very broad Ca-0 feature in the red. The two 'emission 
peaks' observed at ~ 4400 and 5200A correspond to the only two regions in the blue that are 
relatively line- free. A similar situation is observed in SN 1998bw (Iwamoto et al. 1998). 

The spectra are characterized by a low temperature, even near maximum, because the rapid 
expansion combined with the relatively low luminosity (from the tail of the light curve we deduce 
that SN 1997ef produced about 0.15M© of 56 Ni, compared to about O.6M in a typical SN la and 
O.7M in SN 1998bw) leads to rapid cooling. Thus the Sill 6355 A line is not very strong. 

Although model CO 100 yields rather good synthetic spectra, it still fails to reproduce the 
observed large width of the OI- Ca II feature in the only near-maximum spectrum that extends 
sufficiently far to the red (5 Dec 1997). An improvement can be obtained by introducing an 
arbitrary flattening of the density profile at the highest velocities. 

Full details of the spectrum synthesis calculations, including insights on the density structure 
and the abundance in the ejecta will be given in a separate paper (Mazzali et al., 1999, in 
preparation). 

5. DISCUSSION 

5.1. Possible Aspherical Effects 

We have shown that the light curve, the photospheric velocities, and the spectra of SN 1997ef 
are better reproduced with the hyper-energetic model CO 100 than with the ordinary SN Ic model 
CO60. However, there remain several features that are still difficult to explain with model CO100. 

(1) The observed velocity of Si II decreases much more rapidly than models predict. It is as 
high as ~ 30,000 km s" 1 at the earliest phase, but it gets as low as ~ 3,000 km s _1 around day 50 
(Fig. [?]). We find that it is difficult to get such a rapid drop of the photospheric velocity not only 
in models CO100 and CO60, but also in other models that can reproduce the light-curve shape 
reasonably well. Models with higher energies and/or smaller masses would be able to reproduce 
the fast evolution of the photospheric velocity, but such models would inevitably produce light 
curves with a narrower peak and a faster tail. 

(2) Obviously, the observed light curve declines slower than model CO100 in the tail part, 
and it is also a bit flatter than the model near the maximum part (Fig. ||). Models with lower 
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energies and/or larger masses are able to give improved fits to both the peak and the tail of the 
light curve. But, then it gets very difficult to reproduce the large photospheric velocities observed 
at early times in SN 1997ef. 

This dilemma might be overcome if we introduce multiple components of the light curve from 
different parts of ejecta moving at different velocities. In fact, the discrepancies may be interpreted 
as a possible sign of asphericity in the ejecta: A part of ejecta moves faster than average to form 
the lines at such high velocities at early phases, while the other part of ejecta expands with a 
lower velocity so that the low-velocity Si II line comes up at later epochs. Having a low-velocity 
component would also make it easier to reproduce the slow tail. 

(3) Extensive mixing of 56 Ni is required to reproduce the short rise time of the light curve. 
According to hydrodynamical simulations of the Rayleigh- Taylor instability in the ejecta of 
envelope-stripped supernovae (Hachisu et al. 1991; Iwamoto et al. 1996), large scale mixing can 
not be expected to occur in massive progenitors, because in the core of such massive stars the 
density gradient is not steep enough around the composition interfaces. One possibility to induce 
such mixing in the velocity space is an asymmetric explosion(e.g., Nagataki, Shimizu, & Sato 

1998) . Higher velocity 56 Ni could reach the ejecta surface so that the effect of radioactive heating 
comes up as early as is required from light curve modeling. 

In order to realize higher densities at low velocity regions without increasing the mass of 
ejecta significantly, it may be necessary that the explosion is somewhat aspherical. If the explosion 
is aspherical, the shock would be stronger and the material would expand at a larger velocity 
in a certain direction, while in its perpendicular direction, on the other hand, the shock would 
be weaker, ejecting lower- velocity material (e.g., Hoflich et al. 1999). The density of the central 
region could be high enough for 7-rays to be trapped even at advanced phases, thus giving rise to a 
slowly declining tail (see Nakamura et al. 1999a for a discussion of SN 1998bw). In the extremely 
asymmetric cases, the material ejection may happen in a jet-like form. A jet could easily bring 
some 56 Ni from the deepest layer to the surface of high velocity. Detailed spectral analysis of 
observed spectra for different epochs are necessary to investigate this issue further. 

5.2. Gamma-Ray Bursts/Supernovae Connection and SN 1997ef 

There have been an increasing number of candidates for the gamma-ray burst(GRB) /supernova 
connection, including GRB980425/SN1998bw (Galama et al. 1998; Iwamoto et al. 1998; Iwamoto 

1999) , GRB970514/SN1997cy (Germany et al. 1999; Turatto et al. 1999), GRB980910/SN1999E 
(Thorsett & Hogg 1999). Two other high-z GRBs may also be associated with supernovae: 
GRB980326 (Bloom et al. 1999) and GRB970228 (Reichart 1999; Galama et al. 1999). The 
optical transients of these GRBs showed significant reddening and temporal slow down (even with 
a second maximum) in their late light curves, which can be fitted by the early power-law decay 
plus the red-shifted light curve of SN 1998bw. 
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As noted in §1, a possible connection between SN 1997ef and GRB971115 has been suggested 
(Wang & Wheeler 1998). Recently another SN Ic, 1998ey, showed a spectrum with very broad 
features, very similar spectra to that of SN1997ef on Dec 17 (Garnavich et al. 1998); but no GRB 
counterpart has been proposed for SN 1998ey. Although this may cast some doubt on the general 
association between hypernovae and GRBs, it must be noted that both SNe 1997ef and 1998ey 
were less energetic events than SN 1998bw. It is possible that a weaker explosion is less efficient 
in collimating the 7-rays to give rise to a detectable GRB (GRB980425 was already quite weak in 
gamma-rays compared to the average GRBs) , or that some degree of inclination of the beam axis 
to the line of sight results in a seemingly weaker supernova and in the non-detection of a GRB. 
Only the accumulation of more data will allow us to address these questions. 

5.3. The Mass of Ejected 56 Ni 

For the study of the chemical evolution of galaxies, it is important to know the mass of 56 Ni, 
M( Ni), synthesized in core-collapse supernovae as a function of the main-sequence mass M ms of 
the progenitor star (e.g., Nakamura et al. 1999b). From our analysis of SN 1997ef, we can add a 
new point on this diagram. 

We evaluate the uncertainty in our estimates of M( 56 Ni) and M ms . We need 0.15 Mq of 56 Ni 
to get a reasonable fit to the light curve of SN 1997ef at a distance D = 52.3 Mpc. The expected 
10% uncertainty in the distance leads to a 20% uncertainty in the 56 Ni mass, i.e., M( 56 Ni) = 0.15± 
0.03 Mq. The distribution of 56 Ni affects the peak luminosity somewhat, but the effect is found 
to be much smaller than that of the uncertainty in the distance. A 10 Mq C+O star corresponds 
to a M ms = 30 - 35 Mq, but the uncertainty involved in the conversion of the core mass to M ms 
may involve a larger uncertainty if the progenitor undergoes a close binary evolution. 

Figure [H] shows M( 56 Ni) against M ms obtained from fitting the optical light curves of SNe 
1987A, 1993J, and 19941 (e.g., Shigeyama & Nomoto 1990; Nomoto et al. 1993, 1994; Shigeyama 
et al. 1994; Iwamoto et al. 1994; Woosley et al. 1994; Young, Baron, & Branch 1995). The 
amount of 56 Ni appears to increase with increasing M ms of the progenitor, except for SN II 1997D 
(Turatto et al. 1998). 

This trend might be explained as follows. Stars with M ms < 25 Mq form a neutron star, 
producing ~ 0.08 ± 0.03 M Q 56 Ni as in SN lib 1993J, SN Ic 19941, and SN 1987A (although 
SN 1987A may be a borderline case between neutron star and black hole formation). Stars with 
M ms > 25 Mq form a black hole (e.g., Ergma & van den Heuvel 1998); whether they become 
hypernovae or ordinary SNe II may depend on the angular momentum in the collapsing core. For 
SN 1997D, because of the large gravitational potential, the explosion energy is so small that most 
of 56 Ni fell back onto a compact star remnant; the fall-back might cause the collapse of the neutron 
star into a black hole. The core of SN II 1997D might not have a large angular momentum, 
because the progenitor had a massive H-rich envelope so that the angular momentum of the core 
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might have been transported to the envelope possibly via a magnetic-field effect. Similarly, a 
negligible amount of ejection of 56 Ni in black hole formation has recently been suggested for X-ray 
Nova Sco (GRO J1655-40), where the companion star of the black hole seems to be enriched with 
S, Si, Mg, and O but not Fe (Israelian et al. 1999). Hypernovae such as SNe 1998bw, 1997ef, and 
1997cy might have rapidly rotating cores owing possibly to the spiraling-in of a companion star in 
a binary system. The outcome certainly depends also on mass-loss rate and binarity. 

As noted in §5.2, it has been claimed that the optical afterglows of GRB's 980326 and 970228 
are better reproduced if a red-shifted light curve of SN 1998bw is superposed on a power-law light 
component (Bloom et al. 1999; Reichart 1999; Galama et al. 1999). A question arising from 
these two examples is whether the supernovae associated with GRBs have a uniform maximum 
luminosity, i.e., whether 0.7 M & 56 Ni production as in SN 1998bw is rather common or not. 
However, the present study of SN 1997ef shows that the 56 Ni mass and thus intrinsic maximum 
brightness of SN 1997ef is smaller than in SN 1998bw by a factor of 4 - 5 (see the next subsection). 
We certainly need more examples for defining the luminosity function and the actual distribution 
of masses of 56 Ni produced in supernovae/hypernovae. 

5.4. Possible Evolutionary Scenarios 

Here we classify possible evolutionary paths leading to C+O star progenitors. In particular, 
we explore the paths to the progenitors that have rapidly rotating cores with a special emphasis, 
because the explosion energy of hypernovae may be extracted from rapidly rotating black 
holes (Blandford & Znajek 1977). 

(1) Case of a single star: If the star is as massive as M ms > 40 Mq, it could lose H and He 
envelopes in a strong stellar wind (e.g., Schaller et al. 1992). This would be a Wolf-Rayet star. 

(2) Case of a close binary system: Suppose we have a close binary system with a large 
mass ratio. In this case, the mass transfer from star 1 to star 2 inevitably takes place in a 
non-conservative way, and the system experiences a common envelope phase where star 2 is 
spiraling into the envelope of star 1. If the spiral- in releases enough energy to remove the common 
envelope, we are left with a bare He star (star 1) and a main-sequence star (star 2), with a reduced 
separation. If the orbital energy is too small to eject the common envelope, the two stars merge 
to form a single star (e.g., van den Heuvel 1994). 

(2-1) For the non-merging case, possible channels from the He stars to the C+O stars are as 
follows (Nomoto, Iwamoto, & Suzuki 1995). 

(a) Small-mass He stars tend to have large radii, so that they can fill their Roche lobes more 
easily and lose most of their He envelope via Roche lobe overflow. 

(b) On the other hand, larger-mass He stars have too small radii to fill their Roche lobes. 
However, such stars have large enough luminosities to drive strong winds to remove most of the 
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He layer (e.g., Woosley, Langer, & Weaver 1995). Such a mass-losing He star would corresponds 
to a Wolf-Rayet star. 

Thus, from the non-merging scenario, we expect two different kinds of SNe Ic, fast and slow, 
depending on the mass of the progenitor. SNe Ic from smaller mass progenitors (channel a) show 
faster light-curve and spectral evolutions, because the ejecta become more quickly transparent to 
both gamma-ray and optical photons. The slow SNe Ic originate from the Wolf-Rayet progenitors 
(channels b and 1). The presence of both slow and fast SNe Ib/Ic has been noted by Clocchiatti 
& Wheeler (1997). 

(2-2) For the merging case, the merged star has a large angular momentum, so that its 
collapsing core must be rotating rapidly. It would lead to the formation of a rapidly rotating black 
hole from which possibly a hyper-energetic jet could emerge. If the merging process is slow enough 
to eject H and He envelopes, the star would become a rapidly rotating C+O star. Such C+O 
stars are the candidates for the progenitors of Type Ic hypernovae like SNe 1997ef and 1998bw. If 
a significant amount of H-rich (or He) envelope remains after merging, the rapidly rotating core 
would lead to a hypernova of Type Iln possibly like SN 1997cy (or Type lb). 

6. CONCLUSIONS 

We have shown that the photospheric velocities and the spectra of SN 1997ef are much 
better reproduced by the hyper-energetic model CO100 than by the ordinary SN Ic model CO60. 
The model parameters determined for CO 100 are Ek = 8 x 10 51 erg, Mqo = WM Q (which 
corresponds to the C+O core of a 30-35 M© star), and M( 56 Ni)= 0.15 M & . The compact star 
remnant of CO100 is as massive as M rcm ~ 2.4 M & , thus possibly being a black hole. This high 
explosion energy would be extracted from the rapidly rotating black hole. 

For SN 1997ef, M co , M( 56 Ni) and E K are all slightly smaller than for SN 1998bw, but 
SN 1997ef can certainly be regarded as a hypernova in terms of the kinetic energy of explosion. 
Therefore, we suggest that SNe 1997ef, 1998ey, and 1998bw form a new class of hyper-energetic 
Type Ic supernovae, which we call hypernovae. They are distinguished by their large kinetic 
energies, 8-60 times larger than in ordinary supernovae. 

The smaller line velocities at advanced phases and the flatter light curve tail of SN 1997ef 
than the models predict may suggest the presence of a low-velocity, relatively dense core, while 
its higher line velocities at early phases imply the presence of a yet-higher- velocity component of 
ejecta. These are very difficult to be reconciled with any spherically symmetric models, even with 
the high-energy spherical model CO 100. This discrepancy between models and the observations, 
as well as the extensive mixing of 56 Ni required to explain the early rise of the light curve, seems 
to indicate that the explosion of SN 1997ef was at least somewhat aspherical. 
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Table 1. Parameters of the CO star models 



model 


C+O core mass (M ) 


ejecta mass (M Q ) 


56 Ni mass (M ) E K (10 51 erg) 


C021 


2.1 


0.9 


0.07 1 


CO60 


6.0 


4.4 


0.15 1 


CO100 


10.0 


7.6 


0.15 8 


C0138 


13.8 


10.8 


0.7 ~ 30 



Table 2. Predicted yields of SN1997ef (M ) 



model 


c 


O 


Si 


S 


Ca 


Fe 


CO60 


5.2xl0~ 2 


3.0 


0.10 


3.7xl0~ 2 


5.7xl0~ 3 


0.16 


CO100 


0.58 


5.6 


0.42 


0.19 


2.5xl0" 2 


0.19 


model 


44 Ti 


56 Ni 


57 M 








CO60 


2.1xl0" 4 


0.15 


5.7xl0~ 3 








CO100 


4.5xl0" 5 


0.15 


5.7xl0" 3 









Table 2. Parameters of the synthetic spectra 



date 


epoch 


L 


Vph 


VSUI 


log pph 


Mass 


T eff 


Tbb 


Brnod 




V ba 


BC 






days 


erg s _1 


km s 


Rq 


g cm ° 


Mq 


K 


K 












29 Nov 


9 


42.17 


15500 


19072 


-12.65 


0.71 


6123 


7666 


17.45 


16.75 


16.7 


0.28 


-16.700 


5 Dec 


15 


42.19 


9500 


13962 


-12.30 


3.02 


6128 


9407 


17.35 


16.63 


16.5 


0.35 


-16.750 


17 Dec 


27 


42.24 


7500 


8011 


-12.67 


4.79 


5291 


6697 


17.70 


16.59 


16.6 


0.26 


-16.875 
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SNe Ic near maximum 
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Fig. 1.- 



Observed spectra of Type Ic supernovae 1998bw, 1997ef, and 19941. 
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Fig. 2. — Absolute magnitudes of Type Ic supernovae: the ordinary SN Ic 19941 (Richmond et al. 
1996a, b), the hypernova SN 1998bw (Galama et al. 1998), and the proposed hypernova SN 1997ef. 
The dashed line indicates the 56 Co decay rate. 
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Fig. 3. — Density distributions against the velocity of homologously expanding ejecta for CO60 
and CO100. 
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Fig. 4. — Chemical composition of model CO60 plotted against the expansion velocity. 
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Fig. 5.- 



Chemical composition of model CO100 plotted against the expansion velocity. 
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Fig. 6. — Calculated Visual light curves of CO60 and CO 100 compared with that of SN 1997ef. 
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Fig. 7. — Evolution of the calculated photospheric velocities of CO60 and CO 100 (solid lines) 
compared with the observed velocities of the Si II 634.7, 637.1 nm line measured in the spectra at 
the absorption core. 
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Fig. 8. — Observed spectra of SN 1997ef (bold lines) and synthetic spectra computed using 
model CO60. The line features seen in the synthetic spectra are much too narrow compared 
with observations. 
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Fig. 9. — Comparison between the observed spectra of SN 1997ef (bold lines) and synthetic 
spectra computed using model CO100 (fully drawn lines). The fits are much improved with CO100 
compared with the ones with CO60. 
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Fig. 10.— Ejected 56 Ni mass versus the main sequence mass of the progenitors of several bright 
supernovae obtained from light curve models. 



